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Abstract 

Models  for  predicting  the  constitutive  behavior  of  spatially-heterogeneous  microstractures  such  as  discontinuously-reinforced  aluminum  (DRA) 
and  other  metallic  matrix  composites  based  on  unit  cell  approaches  generally  do  not  incorporate  higher-order  microstructural  features  such  as  degree 
of  homogeneity  and  spatial  anisotropy  of  the  reinforcement  phase.  Moreover,  more  complex  numerical  models  rarely  encompass  the  volumes  of 
material  necessary  to  ensure  statistical  relevance.  The  present  contribution  offers  an  alternative  approach  for  quantifying  and  then  incorporating  the 
microstructural  homogeneity  of  these  materials  within  an  elastic-plastic  finite  element  code.  An  attempt  is  made  to  model  both  the  micromechanical 
length  scale  associated  with  the  individual  reinforcement  particles  and  the  microstructural  length  scale  associated  with  their  spatial  distribution, 
at  a  greatly-reduced  computational  expense,  by  using  a  volume- averaged,  discretized  approach.  A  key  assumption  in  this  approach  is  that  below 
the  length  scale  of  the  discretization,  the  micro  structure  can  be  modeled  by  a  uniform  array  of  reinforcement  particles.  The  effect  of  the  level 
of  discretization  on  predictions  of  microstructure -property  relationships  are  not  addressed  directly  in  the  present  work.  Detailed  comparisons  of 
the  present  model  with  discretely-modeled  matrix-particle  geometries  will  therefore  form  the  basis  of  a  subsequent  publication.  Nevertheless, 
several  micro  structure-property  relationships  are  developed  which  reveal  empirical  relationships  between  microstructural  homogeneity  and  the 
elastic-plastic  response. 

Published  by  Elsevier  B.V. 
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1.  Introduction 

Microstructural  details  such  as  degree  of  homogeneity 
and  spatial  anisotropy  of  the  reinforcement  phase  can  have 
a  profound  effect  on  the  mechanical  properties  of  spatially- 
heterogeneous  microstructures  such  as  discontinuously- 
reinforced  aluminum  (DRA)  and  other  metallic  matrix  compos¬ 
ites  [1-6].  A  number  of  previous  attempts  have  been  made  to 
incorporate  different  levels  of  microstructural  detail  into  mate¬ 
rial  property  models  for  discontinuously-reinforced  composite 
materials  [2,7-21].  As  computational  power  has  continued  to 
increase,  numerically-based  2D  and  3D  finite  element  methods 
have  become  increasingly  common.  Over  the  last  15  years, 
there  has  also  been  a  progression  in  size  and  complexity  of  the 
models;  from  the  early  periodic  unit  cell  models  of  Christman 
et  al.  [17]  and  Tvergaard  [21]  through  to  more  recent  work 
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on  2D  [11,22-25]  and  3D  [15,18,22,26,27]  elastic-plastic 
analyses  using  multi-particle  models,  including  the  most  recent 
microstructure-based  approaches  [8,9,11,18,23-25].  For  exam¬ 
ple,  Chawla  et  al.  [18]  used  manual  serial-sectioning  to  produce 
a  microstructurally  accurate  3D  model  of  a  2080/SiC/20p 
aluminum  composite  containing  14  SiC  particles.  In  order  to 
faithfully  reproduce  the  individual  particle  morphologies,  a  total 
of  76,000  elements  was  required.  Other  multi-particle  models 
typically  contain  between  10  and  100  particles  and  require 
around  104-105  elements.  If  this  is  extrapolated  to  the  numbers 
of  particles  needed  to  statistically  represent  microstructural 
heterogeneities  that  operate  on  much  larger  length  scales,  e.g. 
particle  clustering,  then  these  “brute  force”  approaches  to  finite 
element  analysis  quickly  approach  the  practical  limit  of  current 
computational  capabilities.  For  example,  microstructural  het¬ 
erogeneities  in  DRA  have  been  reported  on  the  scale  of  several 
millimeters  [1,28].  In  order  to  model  the  105-106  particles 
contained  in  a  1  mm  cube  of  such  a  material,  a  brute  force 
approach  would  require  108-109  elements,  clearly  beyond  the 
limit  of  current  computational  capabilities. 


0921-5093/$  -  see  front  matter.  Published  by  Elsevier  B.V. 
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The  present  work  therefore  represents  a  departure  from  the 
previous  approaches,  in  that  it  attempts  to  reduce  the  com¬ 
plexity  and  size  of  the  overall  model  by  developing  a  generic 
framework  both  to  quantify  the  spatial  homogeneity  and  then 
to  directly  model  the  influence  of  spatial  homogeneity  on 
mechanical  properties,  using  a  finite  element  scheme  coupled 
with  micromechanical  (analytical)  modeling.  The  approach 
is  computationally  inexpensive,  and  therefore  is  capable  of 
modeling  significant  amounts  of  material,  in  either  2D  or  3D. 
Microstructural  details  such  as  clustering  and/or  homogeneity 
can  be  incorporated  into  this  modeling  scheme  without  resorting 
to  unrealistic  approximations  such  as  periodic  boundary  condi¬ 
tions  [15,26]  or  artificial  symmetry  bounds  [11,27].  The  ability 
to  quantify  and  incorporate  statistical  variations  in  reinforce¬ 
ment  volume  fraction  (i.e.  spatial  homogeneity)  into  models  for 
mechanical  behavior  is  seen  as  a  significant  step  forward  in  the 
ongoing  effort  to  understand  the  effects  of  spatial  homogeneity 
on  mechanical  performance  of  extruded  DRA  materials. 

2.  Quantifying  spatial  homogeneity  using  the  MSAAF 
technique 

The  present  work  concentrates  on  predicting  the  elastic  and 
small-strain  plastic  response  of  DRA  microstructures  that  dis¬ 
play  varying  levels  of  spatial  homogeneity,  due  to  differences 
in  their  processing  histories.  The  Multi-Scalar  Analysis  of  Area 
Fractions  (MSAAF)  technique,  originally  developed  by  Spowart 
et  al.  [29],  is  used  to  quantify  the  spatial  homogeneity  of  the 
different  DRA  microstructures  by  dividing  up  the  microstruc¬ 
ture  into  local  grid  squares  and  then  measuring  the  coeffi¬ 
cient  of  variation  of  local  area  fractions  on  a  sectioning  plane, 
CV(Af)  =  ctAf/Af ,  where  OAf  is  the  standard  deviation  of  the 
individual  local  area  fractions  on  the  plane  and  Af  is  their  statis¬ 
tical  mean.  For  a  microstructure  comprising  a  Poisson  random 
distribution  of  mono-sized  circular  particles  of  diameter  dp,  the 
statistical  parameter  should  vary  with  grid  size,  Q ,  according  to 
Eq.  (1)  for 


For  Q<dv ,  however,  it  is  found  that  CV(Af)  depends  only  on  the 
mean  area  fraction  of  the  original  image,  regardless  of  the  spatial 
distribution  of  particles.  In  this  case,  Eq.  (2)  is  more  appropriate 

/ 1  _  Af\0,5 

CV(Af)|lim  e/dp^0  =  (2) 

Fig.  1  shows  MSAAF  data  obtained  from  a  model  microstruc¬ 
ture  generated  using  the  RSA  algorithm  [29]  comprising  10,000 
randomly-distributed  non- overlapping  mono-sized  disks  with  an 
overall  area  fraction  of  Af  =  0.022.  There  is  excellent  agreement 
between  the  calculated  values  of  CV(Af)  and  the  predictions 
of  Eqs.  (1)  and  (2)  for  an  equivalent  Poisson  random  distri¬ 
bution  of  points.  The  slope  of  the  straight-line  portion  on  this 
log-log  plot  is  close  to  —  1 .0,  due  to  the  random  nature  of  the 
micro  structure.  The  technique  provides  microstructural  infor¬ 
mation  over  many  length  scales,  from  roughly  the  particle  size 


Fig.  1.  MSAAF  results  showing  variation  of  CV(Af)  with  normalized  length 
scale,  Q/dp  for  random,  non-overlapping  mono-sized  disks  with  overall 
Af  =  0.022.  Also  plotted  are  lines  representing  Eqs.  (1)  and  (2)  for  the  equiv¬ 
alent  Poisson  random  distribution.  The  small  black  triangle  on  the  lower  right 
horizontal  axis  represents  the  1%  homogeneous  length  scale,  =  580  dp. 

upwards,  however,  it  is  often  convenient  to  characterize  the  level 
of  spatial  heterogeneity  in  the  microstructure  by  assigning  a  sin¬ 
gle  scalar  parameter,  the  Homogeneous  Length  Scale  or  Lh  .  This 
parameter  (with  units  of  length)  represents  the  microstructural 
length  scale  beyond  which  the  micro  structure  is  homogeneous 
to  within  a  prescribed  range.  For  example,  one  may  define  the 
parameter  L^01  as  the  length  scale  beyond  which  the  local  sta¬ 
tistical  variation  in  area  fraction  is  no  greater  than  1%,  indicated 
by  the  triangular  symbol  in  Fig.  1.  For  this  micro  structure  of 
non-overlapping  circular  discs  with  an  overall  area  fraction  of 
0.022,  L^01  =  580  x  dp,  where  dv  is  the  diameter  of  the  (mono¬ 
sized)  particles,  however,  for  experimental  micro  structures  L^01 
is  typically  measured  in  jim  or  mm.  When  comparing  differ¬ 
ent  microstructures  with  similar  area  fractions,  the  smaller  the 
homogeneous  length  scale,  the  more  homogeneous  (or  less  clus¬ 
tered)  the  microstructure.  Alternatively,  one  may  consider  the 
slope,  s  as  an  indicator  of  microstructural  heterogeneity.  Steeper 
slopes  (i.e.  more  negative)  suggest  more  uniform  microstruc¬ 
tures,  whereas  shallower  slopes  (i.e.  less  negative)  suggest  more 
clustered  or  heterogeneous  microstructures.  Generally,  this  is  a 
less  sensitive  descriptor  of  homogeneity  than  Lh- 

2.7.  Extension  of  MSAAF  for  anisotropic  micro  structures 

It  was  anticipated  that  extrusion  of  the  DRA  would  lead  to 
microstructural  anisotropy  in  the  longitudinal  plane.  In  order 
to  include  microstructural  anisotropy  in  the  measurement  of 
spatial  homogeneity,  the  original  MSAAF  technique  was  mod¬ 
ified  to  allow  for  anisotropic  resampling.  Instead  of  using 
isotropic  resampling,  whereby  pixels  in  the  new  image  replace 
square  arrays  of  pixels  in  the  original  image  (i.e.  1  x  1,  2  x  2, 
4x4,  etc.)  anisotropic  resampling  is  used,  where  single  pix¬ 
els  replace  strips  of  pixels  in  the  original  image,  i.e.  1x2, 
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1x4,  1  x  8,  etc.  for  strips  oriented  horizontally  on  the  image 
(parallel  to  the  extrusion  direction),  and  2x1,  4x1,  8x1, 
etc.  for  strips  oriented  vertically  (transverse  to  the  extrusion 
direction). 

The  anisotropic  analysis  using  oriented  pixel  strips  enables 
differentiation  between  the  two  orthogonal  directions  in  the 
image.  The  means  and  standard  deviations  of  all  the  gray- 
levels  in  each  strip  are  calculated  as  for  the  isotropic  case, 
but  are  now  used  to  compute  two  separate  coefficients  of  vari¬ 
ation,  CVi  for  horizontal  strips  and  CV2  for  vertical  strips. 
This  allows  two  different  MSAAF  curves  to  be  plotted,  and 
leads  to  the  definition  of  two  separate  (in-plane)  homoge¬ 
neous  length  scale  parameters,  Lhi  and  Lh2-  For  randomly- 
distributed  non-overlapping  mono- sized  disks,  the  slopes  of 
the  anisotropic  MSAAF  curves  will  be  close  to  —0.5  (cf. 
—  1.0  for  the  isotropic  case)  which  makes  extrapolation  of 
the  straight-line  portion  of  the  curve  to  a  CV  of  0.01  some¬ 
what  impractical.  For  this  reason,  Lh  is  typically  defined  at 
a  CV  of  10%  (i.e.  L^jandL^)  for  anisotropic  microstruc¬ 
tures  such  as  these.  In  addition,  the  ratio  of  the  two  homoge¬ 
neous  length  scales  is  a  measure  of  the  amount  of  anisotropy 
in  the  micro  structure,  and  is  therefore  another  microstructural 
descriptor  which  may  be  related  to  processing  conditions  (see 
Section  3). 


3.  Experimental 

Spherical,  gas-atomized  matrix  alloy  powders  of  nominal 
composition  Al-0.27  Cu-0.26  Fe-0.97  Mg-0.56  Si  (wt.%,  6061 
aluminum  alloy)  were  classified  into  three  powder  lots  with  aver¬ 
age  particle  sizes  of  26.4,  42.0  and  108.6  p,m.  These  powders 
were  then  blended  with  high  purity  F600  grade  SiC  (average  par¬ 
ticle  size  13.4  |jim)  to  produce  three  different  composite  materi¬ 
als  with  increasing  levels  of  spatial  heterogeneity,  as  the  particle 
size  ratio  (PSR)  is  increased  from  2:1  to  3:1  to  8:1  [1].  Fig.  2 
shows  optical  micrographs  of  longitudinal  sections  through  typ¬ 
ical  composite  microstructures  for  each  of  the  three  different 
materials. 

Five  separate  regions  were  carefully  selected  from  each 
longitudinally- sectioned  DRA  microstructure,  in  order  to  pro¬ 
vide  adequate  statistics  for  the  measurement  of  spatial  homo¬ 
geneity,  and  also  to  ensure  that  several  duplicate  images  with 
the  same  overall  area  fraction  were  available  for  numerical  mod¬ 
eling.  The  area  of  each  image  was  667  |mm  x  667  p,m,  and  the 
mean  area  fraction  of  the  images  was  measured  as  0.275  zb  0.003. 
Additionally,  particle  aspect  ratios  in  the  longitudinal  section 
plane  were  measured  using  quantitative  image  analysis  software 
(Image  Processing  Toolkit,  version  3.0)  [30].  The  average  aspect 
ratio  of  the  F-600  grade  SiC  particles,  s,  was  calculated  to  be  2.0. 


Fig.  2.  Optical  micrographs  showing  longitudinal  sections  through  extruded  DRA  composites:  (a)  PSR  =  2:1,  (b)  PSR  =  3:1,  (c)  PSR  =  8:1.  Extrusion  direction  is 
horizontal.  Note  decreasing  spatial  homogeneity  and  increasing  microstructural  anisotropy  with  increasing  PSR. 
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4.  Modeling  of  composite  mechanical  response 

A  combined  model  was  developed  to  describe  the  effects  of 
microstructure  on  the  composite  mechanical  response  for  the  dif¬ 
ferent  DRA  materials.  Composite  behavior  at  the  length  scale 
of  the  reinforcement  particles  was  modeled  using  a  microme¬ 
chanical  approach  based  on  the  Eshelby  equivalent  inclusion 
method  [31,32]  for  predicting  the  elastic  and  small-strain  plas¬ 
tic  response  of  idealized  two-phase  materials  containing  ellip¬ 
soidal  inclusions.  The  overall  composite  behavior  was  modeled 
using  a  finite  element  code  to  discretize  the  microstructure  into 
individual  domains  which  had  constitutive  properties  given  by 
the  Eshelby  model.  The  combined  model  thus  contains  both 
micromechanical  details  at  shorter  length  scales  and  microstruc- 
tural  information  at  longer  length  scales,  thereby  providing  a 
fairly  efficient  framework  for  modeling  the  different  spatial 
distributions  of  the  reinforcement  particles  without  requiring 
the  level  of  complexity  associated  with  a  brute  force  meshing 
approach. 

4.1.  Micromechanical  modeling  using  the  Eshelby 
equivalent  inclusion  method 

The  Eshelby  equivalent  inclusion  method  is  a  technique 
whereby  the  elastic  behavior  (and  to  a  limited  extent,  the  plastic 
behavior)  of  a  two-phase  material  containing  ellipsoidal  inclu¬ 
sions  can  be  predicted,  based  on  the  calculation  of  internal 
stresses  that  arise  due  to  a  generalized  misfit  strain  between 
the  shapes  of  the  two  constituent  phases,  matrix  and  inclusion.1 
For  example,  during  elastic  loading,  the  misfit  may  arise  due 
to  a  stiffer  inclusion  wishing  to  deform  to  a  lesser  extent  than 
a  more  compliant  matrix  material  under  the  same  applied  load, 
whereas  during  plastic  straining,  the  misfit  might  arise  due  to 
plastic  deformation  of  the  matrix  whilst  the  inclusion  remains 
elastic.  The  Eshelby  analysis  shows  that  for  a  generalized  misfit 
strain,  £*T,  the  internal  stress  within  an  ellipsoidal  inclusion,  cr\, 
is  given  by  [33] 

cri  =  Cm(S  -  /){(Ci  -  Cm)S  +  Cm}~lCis*T  (3) 

Here  Cm  is  the  stiffness  tensor  of  the  matrix  material,  S  the 
Eshelby  tensor,  I  the  identity  matrix,  and  C\  is  the  stiffness  tensor 
of  the  inclusion.  In  order  to  extend  the  Eshelby  model  to  cope 
with  non-dilute  systems  containing  an  inclusion  volume  fraction 
fi  the  mean-field  approximations  of  Mori  and  Tanaka  [34]  and 
Pedersen  [35]  are  used,  whereby  the  volume-averaged  stresses 
in  each  constituent  phase,  (cr)m  and  (cr)i  can  be  given  in  terms 
of  the  equivalent  transformation  strain,  sT  as 

<tf)m  =  ~fCm(S  -  /)eT  (4) 

and 

<tf>i  =  (1  -  f)Cm(S  -  I)eT  (5) 


1  In  what  follows,  the  particles  are  referred  to  as  inclusions,  in  keeping  with 
Eshelby’s  original  treatment. 


The  composite  elastic  stiffness  tensor,  Cc,  can  then  be  obtained 
directly  from  the  constituent  volume- averaged  stresses,  (cr)m 
and  (a) i  as 

Cc  =  [C"1  -  /{(Ci  -  Cm)[S  -  f(S  -  /)]  -  Cm}-1 

(Cl  -  C^C-'f 1  (6) 

A  freely-available  computer  code  [36]  was  used  to  evaluate 
this  expression  for  Cc  at  volume  fractions  between  0  (pure 
matrix)  and  1  (pure  reinforcement)  in  order  to  provide  mate¬ 
rial  property  data  for  input  into  the  finite  element  model  as 
engineering  constants.  Due  to  symmetry,  the  stiffness  ten¬ 
sor  for  transversely-isotropic  materials  such  as  the  extruded 
DRA  composite  reduces  to  seven  independent  components. 
Expressed  in  terms  of  engineering  constants,  these  components 
are En,E22  =  E33,  v12  =  v13,  v2i  =  v31,  v23=v32 ,  G\2  =  G\3  and 
G23  =  G32 ,  where  the  one-direction  is  chosen  to  be  parallel  to  the 
extrusion  axis,  and  the  (equivalent)  two-  and  three-directions  are 
transverse  to  the  extrusion  axis. 

4.1.1.  Composite  yield  criterion 

The  Tresca  or  maximum  shear  stress  yield  criterion  used  for 
the  composite  can  be  generally  stated  as  Aai2  =  2k ,  where  Aai2 
is  the  maximum  deviatoric  stress  measured  in  the  (1-2)  plane, 
and  k  is  the  yield  stress  measured  in  pure  shear.  For  uniaxial 
tension  in  the  axial  direction  under  an  applied  load  <7A,  o\  -  crA 
and  a i  =  o2  =  0,  therefore 

A  cri2  =  crA  =  cry  =  2k  (7) 

For  a  transversely-isotropic  material  where  the  two-  and  three- 
axes  are  equivalent 

cry  =  (dim  -  02m)  =  (film  ~  &3m)  =  A<7]2m  (8) 

Here  crym  is  the  matrix  uniaxial  yield  stress  and  dim,  <72m, 
and  <73m  are  given  by  the  sum  of  the  mean  internal  misfit 
stresses  arising  due  to  the  stiffness  mismatch  between  parti¬ 
cle  and  matrix,  (cri),  and  the  appropriate  component  of  the 
applied  stress,  aA.  For  uniaxial  tension  in  the  axial  direction, 
crA  >  0,  aA  =  oA  =  0,  therefore  the  composite  yield  criterion 
based  on  volume-averaged  matrix  stresses  becomes 

tfYlc  =  m  ((cri)  m  (°r2>m)  (9) 

The  yield  stress  of  the  composite  is  higher  than  that  of  the  matrix 
material  since  (ori)m<  (<72  )m  for  all  particle  aspect  ratios  >1. 
Furthermore,  the  matrix  deviatoric  stress  (Acri2m)  will  scale 
linearly  [33]  with  the  applied  stress,  i.e.  Acri2m  =  P\crA.  The 
scalar  parameter,  P\ ,  effectively  describes  the  difference  in  load 
partitioning  between  the  axial  and  transverse  directions  in  the 
composite,  averaged  over  all  the  matrix  volume.  Accordingly, 
the  yield  criterion  for  the  matrix  (and  hence  the  composite)  can 
be  written 


For  the  particle-reinforced  composite  system  Al/SiC/27.5p,  P\ 
evaluates  to  0.724  for  axial  loading  of  inclusions  with  aspect 
ratio  s  =  2,  increasing  the  yield  stress  in  the  axial  direction  by 
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around  35%,  compared  to  the  un-reinforced  matrix.  Similar  cal¬ 
culations  can  also  be  performed  for  transverse  loading,  however, 
due  to  the  aspect  ratios  of  the  particles,  load-sharing  in  the  trans¬ 
verse  direction  is  slightly  less  efficient  than  in  the  axial  direction, 
therefore  aY2C  =  aY3C<aYic.  For  the  purposes  of  the  present 
work,  the  yield  surface  is  assumed  to  be  an  ellipsoid,  with  major 
axes  given  by  aY ic,  <?Y2c  and  aY3c- 


4.1.2.  Composite  strain-hardening  behavior 

This  treatment  is  also  applicable  when  the  misfit  strain  arises 
as  a  result  of  matrix  plasticity.  With  the  assumption  of  perfectly- 
plastic  deformation  at  constant  volume  and  a  uniaxial  plastic 
strain  of  sfm,  the  transformation  strain  for  an  isolated  inclusion 
is  given  by 


s*T  = 


£lm 

iopl 

2fclm 

icpl 

2£lm 


(ID 


This  provides  another  matrix  deviatoric  stress  term,  Aof2 m? 
which  depends  linearly  on  the  level  of  matrix  plastic  strain. 
Using  the  approximation  sf  =  (1  —  f)£m  [33],  it  follows  that 
A  erf 2m  =  —  Wg,  where  sf  is  the  average  plastic  strain  along  the 
axial  direction,  and  W  has  units  of  MPa/p,s  (for  Al/SiC/27.5p, 
W=  17.7  GPa).  With  a  perfectly-plastic  matrix,  linear  strain¬ 
hardening  is  predicted  for  the  composite,  i.e. 


<rYm  =  -Wsf  +  PaA 


(12) 


In  this  approximation,  the  only  contribution  to  strain-hardening 
is  from  load  partitioning  between  the  elastic  particles  and  the 
perfectly  plastic  matrix.  However,  any  hardening  (or  soften¬ 
ing)  that  may  occur  in  the  matrix  during  deformation  will  have 
additional  effects.  For  example,  parabolic  work-hardening  as 
often  seen  in  aluminum  alloys  would  dramatically  increase 
the  strain-hardening  rate.  Conversely,  relaxation  mechanisms 
such  as  interfacial  debonding,  cavitation  [37],  dislocation  move¬ 
ment  [38]  and/or  interfacial  diffusion  [39]  can  lower  the  overall 
strain-hardening  rate.  In  an  attempt  to  better  predict  the  ini¬ 
tial  strain-hardening  response  of  the  composite,  an  additional 
(matrix)  strain-hardening  term  which  accounts  for  the  initial 
work-hardening  of  the  6061  aluminum  alloy  matrix  is  included 
in  Eq.  (12),  to  give  an  overall  strain-hardening  rate  of 


da  _  W  da 

dsf~T  +  d^ 


(13) 


Fig.  3  shows  the  predicted  elastic-plastic  responses  for 
three  different  reinforcement  volume  fractions  of  aligned 
(5  =  2.0)  particles,  using  the  following  handbook  data  [33,40] 
for  the  matrix  and  reinforcement  constituent  properties: 
Ea\  =  70  GPa,  £siC  —  400  GPa,  vai  —  0.33,  vsic  —  0.20,  ay  (6061- 
T6)  =  0.400  GPa,  Aa/Ast  (6061-T6)  =  8.0  GPa. 

Due  the  alignment  of  the  particles,  axial  modulus  and  yield 
strength  is  predicted  to  be  higher  than  transverse  modulus  and 
yield  strength.  Both  properties  increase  with  increasing  rein¬ 
forcement  volume  fraction  due  to  enhanced  load- sharing.  For 


Fig.  3.  Composite  axial  and  transverse  tensile  responses,  as  predicted  by  Eshelby 
model  for  reinforcement  volume  fractions  of  0.1  and  0.6.  Also  shown  are  pre¬ 
dicted  stress-strain  responses  in  the  axial  and  transverse  loading  directions  for 
a  composite  with  overall  volume  fraction  of  0.275,  corresponding  to  the  (hypo¬ 
thetical)  homogeneous  composite  material,  Al-6061 /SiC/27. 5p. 

volume  fractions  less  than  0.5,  the  initial  strain-hardening  rate, 
da/dsVc  increases  with  increasing  volume  fraction  as  load¬ 
sharing  becomes  more  effective  whilst  for  volume  fractions 
greater  than  0.5,  da/d£cl decreases  due  to  the  decreasing  volume 
fraction  of  matrix.  The  Eshelby  model  therefore  predicts  the 
maximum  strain-hardening  response  at  Vf  =  0.5.  The  curve  for 
Vf  =  0.275  corresponds  to  the  Eshelby  prediction  for  a  hypothet¬ 
ical  “homogeneous”  composite  with  the  same  overall  volume 
fraction  as  the  average  area  fraction  measured  in  the  actual 
composite  microstructures,  Af  =  0.275.  Although  it  is  not  a  rig¬ 
orous  upper  bound,  the  “homogeneous”  material  represents  a 
practical  upper  limit  for  composite  mechanical  properties,  since 
microstructural  heterogeneity  has  a  negative  effect  on  mechan¬ 
ical  properties  in  these  materials. 

4.2.  Microstructural  modeling  using  finite  elements 

Since  the  Eshelby  approach  provides  no  spatial  informa¬ 
tion  for  modeling,  a  finite  element  scheme  was  developed  to 
combine  the  spatial  information  contained  in  the  microstruc¬ 
tural  images  with  the  constitutive  responses  of  the  individual 
elements  predicted  by  Eshelby.  Each  element  is  modeled  as  a 
microstructural  domain  consisting  of  a  uniform  spatial  distribu¬ 
tion  of  ellipsoidal  reinforcement  particles,  at  a  particular  volume 
(area)  fraction  calculated  from  the  individual  pixel  gray-levels  in 
the  corresponding  digital  micrograph.  A  custom  pre-processor 
was  developed  in  order  to  correctly  assign  the  256  different 
material  property  definitions  used  in  the  finite  element  code, 
corresponding  to  the  256  available  gray-levels  in  the  image. 
Although  this  is  an  approximation  of  the  true  microstructural 
details  at  the  micromechanical  length  scale,  the  assumption  of 
a  uniform  distribution  of  reinforcement  inside  each  element 
enables  the  model  to  extend  over  a  large  area  of  microstruc¬ 
ture  without  requiring  a  prohibitively  large  number  of  elements. 
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Each  of  the  models  presented  in  this  work  covers  an  area 
of  667  fxm  x  667  p,m,  corresponding  to  16,384  microstructural 
domains  (finite  elements)  in  a  128  x  128  square  mesh.  The  level 
of  discretization  was  chosen  for  computational  expediency,  with 
the  size  of  each  element  being  much  smaller  than  the  homoge¬ 
neous  length  scale.  Each  667  juim  x  667  p,m  model  is  estimated 
to  contain  at  least  1000  particles,  by  stereological  analysis.  It 
was  anticipated  that  there  would  be  significant  amounts  of  scat¬ 
ter  in  the  predicted  stress-strain  curves  from  the  FEM,  owing  to 
the  microstructural  variability  inherent  in  the  composite  mate¬ 
rials.  A  full  3D  solution  was  therefore  considered  (rather  than  a 
2D  plane-stress,  plane-strain  or  axisymmetric  solution)  in  order 
to  provide  more  precise  results.  The  3D  solution  can  also  be 
easily  extended  to  multiple  layers,  for  example  to  model  3D 
microstructural  data  obtained  via  automated  serial- sectioning 
or  computed  X-ray  tomography  techniques,  as  these  datasets 
become  more  readily  available  [41] .  Eight-noded  linear  3D  brick 
elements  (C3D8)  with  equal  side  lengths  were  chosen  through¬ 
out,  with  no  variation  of  properties  with  depth.  During  axial 
and  transverse  loading,  one  side  of  the  mesh  is  subjected  to  a 
uniform  displacement  in  the  loading  direction,  with  the  oppo¬ 
site  side  of  the  mesh  constrained  for  zero  displacement  along 
the  same  direction.  In  addition,  one  corner  of  this  constrained 
side  is  further  constrained  in  all  degrees  of  freedom  to  prevent 
rigid-body  motion  of  the  mesh.  The  analysis  was  carried  out 
using  ABAQUS/Standard  (version  6.4-3)  on  a  single-processor 
machine,  with  total  run-times  less  than  1  h,  providing  elastic- 
plastic  solutions  for  10  separate  strain  steps  up  to  a  total  strain 
of  1.0%. 


4.2.1.  Validation  of  FEM  method 

In  order  to  validate  the  FEM  method  outlined  above,  an 
attempt  was  made  to  compare  output  from  the  model  with 
previously  published  [23,42]  experimental  stress-strain  data 
for  an  extruded  Al-2080/SiC/30p-T6  DRA  material,  of  sim¬ 
ilar  microstructure  to  the  present  DRA.  Optical  micrographs 
were  obtained  from  the  authors,  showing  longitudinal  sections 
through  the  material.  A  507  |mm  x  347  jim  region  of  the  DRA 
microstructure,  containing  the  axial  and  transverse  directions 
was  selected,  digitized,  binarized  (thresholded)  and  resampled 
to  provide  a  95  x  65  rectangular  grid  of  square  finite  elements. 
As  before,  each  element  was  given  properties  according  to 
its  volume  fraction  (gray  level)  following  the  Eshelby  treat¬ 
ment.  The  matrix  and  reinforcement  constituent  properties  used 
for  the  model  were  based  on  Chawla’s  original  measurements: 
£ai  —  74  GPa,  £sic  —  410  GPa,  vai  =  0.33,  vsic  =  0.19,  cry  (2080- 
T6)  =  0.390  GPa,  Aa/Ast  (2080-T6)  =10.8  GPa,  s  =  2.00. 

Fig.  4  shows  there  is  fairly  good  agreement  between  the 
experimental  stress-strain  behavior  of  the  Al-2080/SiC/30p 
material  and  the  prediction  from  the  combined  Eshelby-FEM 
technique.  The  elastic  portion  of  the  curve  is  very  well  matched, 
and  the  strain-hardening  rate  is  similar  over  quite  a  broad  range 
of  total  strain,  up  to  ~1.2%.  Clearly,  as  strain  increases  beyond 
this,  the  assumption  of  linear  work-hardening  in  the  matrix 
becomes  less  appropriate,  since  localized  damage  may  be  occur¬ 
ring  in  the  experimental  material  which  is  not  addressed  in  the 


Fig.  4.  Comparison  between  experimentally-generated  [42]  stress-strain  curve 
for  extruded  Al-2080/SiC/30p-T6  DRA  material  and  results  from  combined 
Eshelby-FEM  (axial  loading). 


present  model.  However,  the  model  appears  appropriate  for  low 
strains  (i.e.  less  than  1%).  There  is  also  a  significant  difference 
in  the  two  curves  in  the  region  around  the  yield  point,  with 
the  yield  strength  of  the  model  being  higher  than  in  the  experi¬ 
ment.  This  is  most  likely  due  to  volume-averaging  of  the  matrix 
stresses  in  the  present  model  which  will  tend  to  overestimate  the 
stress  necessary  to  cause  yielding  [33] .  Additional  complications 
may  arise  due  to  thermal  residual  stresses  from  processing  the 
DRA.  Although  micro-yielding  phenomena  are  not  explicitly 
addressed  in  the  present  model,  their  effect  and  the  associated 
influence  of  thermal  residual  stresses  will  be  investigated  in  a 
future  publication. 

5.  Results 

5.1.  Quantification  of  spatial  homogeneity  in  anisotropic 
micro  structures 


The  anisotropic  MSAAF  procedure  was  used  to  analyze  a 
total  of  15  separate  digital  images,  corresponding  to  five  dif¬ 
ferent  667  p,m  x  667  jxm  regions  of  interest  (ROI’s)  in  each  of 
the  three  different  extruded  DRA  microstructures.  Fig.  5  shows 
the  axial  and  transverse  variations  of  CV(Af)  for  each  material, 
as  a  function  of  strip  length,  Q ,  averaged  over  the  five  different 
ROI’s.  The  plots  also  show  least- squares  curve  fits  to  the  data, 
using  Eq.  (14),  where  a  is  a  geometric  factor  related  to  the  mean 
particle  size  (dv)  and  the  mean  area  fraction  (Af)  and  §  is  the 
slope  of  the  straight-line  portion  on  the  log-log  plot 


CV  ( Af)  = 


Af 


1  -  Af 


+  oiQ 


- 2 £ 


-i  -0.5 


(14) 


The  anisotropic  homogeneous  length  scales,  and  are 
plotted  as  inverted  triangles  on  the  horizontal  axis.  These  are 
obtained  by  extrapolating  the  least-squares  curve  fits  from  Eq. 
(14)  to  CV(Af)  =  0.1.  There  is  a  steady  decrease  in  spatial  homo- 
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(a)  strip  length,  Q  (microns) 


Fig.  5.  Anisotropic  MSAAF  results  for  extruded  DRA  microstructures  in:  (a) 
axial,  and  (b)  transverse  directions.  The  inverted  triangles  on  the  horizontal  axis 
are  the  in-plane  homogeneous  length  scales,  andL^J-  Not  only  is  there  a 
decrease  in  spatial  homogeneity  with  increasing  PSR,  the  spatial  homogene¬ 
ity  for  each  DRA  microstructure  is  lower  in  the  extrusion  direction  than  in  the 
transverse  direction  due  to  microstructural  anisotropy  introduced  by  the  extru¬ 
sion  process.  In  addition,  the  slopes  of  the  linear  portions  on  the  plots  rank 
approximately  with  increasing  PSR,  indicating  decreasing  spatial  homogeneity. 

geneity  with  increasing  PSR,  as  characterized  by  increases  in 
both  anisotropic  homogeneous  length  scales,  and  •  Fur¬ 
thermore,  at  each  PSR,  the  spatial  homogeneity  is  lower  in 
the  extrusion  direction  than  in  the  transverse  direction  due  to 
microstructural  anisotropy  introduced  by  the  extrusion  process. 
The  ratio  Lhi/£h2  can  be  used  as  an  indication  of  the  amount 
of  anisotropy  in  the  microstructure,  and  is  typically  found  to 
increase  with  increasing  PSR. 

Additional  evidence  for  decreasing  spatial  homogeneity  with 
increasing  PSR  is  provided  by  the  slopes  of  the  linear  portions 
on  the  plots  ranking  approximately  with  increasing  PSR.  This  is 
easily  understood,  since  over  all  length  scales  there  will  be  higher 
values  of  CV(Af)  as  the  microstructural  variability  increases. 
However,  the  measured  slopes  are  generally  a  less-sensitive  mea¬ 


Table  1 


Anisotropic  microstructural  descriptors  for  each  extruded  DRA  microstructure 


Particle  size  ratio 

L°h;  (pun) 

£i 

^H2  <Am) 

£2 

2:1 

1250 

-0.545 

950 

-0.558 

3:1 

1590 

-0.535 

1250 

-0.531 

8:1 

6650 

-0.426 

2470 

-0.495 

sure  of  spatial  homogeneity  than  the  change  in  homogeneous 
length  scale,  and  they  are  more  prone  to  experimental  scatter. 
The  averaged  values  for  Lhi  ,  Lh2,  and  §  are  presented  in  Table  1 . 

A  large  number  of  different  microstructural  descriptors  for 
homogeneity  have  been  previously  presented  in  the  literature. 
These  include  descriptors  based  on  Dirichlet  and  other  tessella¬ 
tions  [3,6,43-46],  inter-particle  spacings  [5,8,47],  radial  distri¬ 
bution  functions  [48],  ft-point  correlation  functions  [28,49,50], 
skeletonization  techniques  [51]  and  measurements  of  local  area 
fractions  [1,29,52-54].  In  order  to  develop  the  phenomeno- 


Fig.  6.  Typical  tensile  stress-strain  curves  for:  (a)  axial  loading,  and  (b)  trans¬ 
verse  loading  of  finite  element  models  of  each  of  the  different  classes  of  extruded 
composites  and  homogeneous  model  microstructures.  Although  the  elastic  load¬ 
ing  portions  are  similar,  the  effect  of  microstructure  is  evident  in  the  plastic 
behavior  of  each  material. 
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logical  models  necessary  for  making  the  first  steps  towards 
understanding  the  complex  relationships  between  microstruc- 
ture  and  mechanical  behavior,  it  is  important  to  establish  which 
microstructural  descriptors  are  most  relevant.  The  homogeneous 
length  scales  £ri  and  £r2  which  result  from  the  anisotropic 
MSAAF  analyses  are  offered  as  potentially-useful  additional 
microstructural  descriptors  which  capture  both  homogeneity  and 
anisotropy  information  in  a  quantitative  way. 

5.2.  Effect  of  spatial  homogeneity  on  composite  tensile 
stress-strain  curves 

Fig.  6  shows  the  axial  and  transverse  tensile  stress-strain 
curves  generated  from  the  finite  element  models,  for  each  of  the 
different  classes  of  extruded  composite.  Each  curve  represents 
tensile  data  plotted  from  the  model  which  gave  results  closest  to 
the  average  response  for  all  the  models  in  that  class.  Also  plot¬ 
ted  are  data  from  homogeneous  composite  models,  containing 
elements  with  only  one  material  definition  The  homogeneous 
models  represent  ideal  microstructures,  with  perfectly  uniform 
packing  of  aligned  (s  =  2.0)  particles  at  a  reinforcement  volume 
fraction  of  0.275. 

These  data  show  that  axial  loading  produces  a  higher  com¬ 
posite  stiffness  and  a  higher  yield  strength  (£i  =  1 15  GPa,  o yi  = 
0.558  GPa)  than  in  the  transverse  orientation  (£2  =  103  GPa, 
°ty2  —  0.5 19  GPa),  as  discussed  in  Section  4.  This  is  primarily 
due  to  the  anisotropic  material  properties  that  ensue  from  using 
non-spherical  (i.e.  s  =  2)  particles  in  the  Eshelby  model,  demon¬ 
strating  that  regardless  of  any  possible  additional  microstructural 
effects  (e.g.  due  to  extrusion),  there  is  still  the  potential  to  pro¬ 
duce  anisotropic  mechanical  behavior  in  such  a  composite.  From 
the  elastic  portion  of  the  curves,  there  is  very  little  evidence  of 
any  effect  of  microstructure  on  the  elastic  modulus  of  these  com¬ 
posites,  since  the  curves  lie  almost  exactly  on  top  of  each-other 
up  to  the  point  of  yielding.  There  is,  however,  a  general  trend 
for  the  post-yield  stress-strain  curves  of  the  more  homogeneous 
material  (i.e.  PSR  =  2:1)  to  lie  above  those  for  the  less  homoge¬ 
neous  materials,  PSR  =  3:1  and  PSR  =  8:1,  in  either  orientation. 
The  FEM  results  suggest  that  decreasing  the  spatial  homogene¬ 
ity  of  the  microstructure  has  a  significant  effect  on  the  predicted 
stress-strain  response,  lowering  both  the  yield  strength  and  the 
subsequent  strain-hardening  rate. 

5.3.  Influence  of  spatial  homogeneity  on  elastic  behavior 

Table  2  provides  a  summary  of  the  FEM  results  for  the  elastic 
moduli  of  the  different  composite  microstructures,  along  with 
the  predictions  from  the  Eshelby  treatment.  Although  there  is 


Table  2 

Summary  of  finite  element  results  for  axial  and  transverse  elastic  moduli 
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Fig.  7.  Finite  element  results  for:  (a)  axial,  and  (b)  transverse  composite  yield 
strengths,  plotted  vs.  the  microstructural  descriptors  and  L ^  •  There  is  a 
strong  effect  of  decreasing  microstructural  homogeneity  leading  to  lower  yield 
strengths  in  each  loading  direction,  with  slightly  less  scatter  in  the  yield  strength 


a  clear  difference  between  the  two  loading  directions  due  to 
differences  in  load-sharing,  there  is  clearly  very  little  influence 
of  microstructural  homogeneity  on  E\ .  There  may  be  a  slight 
decrease  in  £2  with  increasing  PSR,  however  this  is  only  a  very 
minor  effect  (i.e.  <2%  over  the  whole  range  of  microstructural 
homogeneities).  This  constancy  of  elastic  modulus  is  expected 
from  the  results  of  previous  numerical  studies  of  whisker- 
and  particle-reinforced  composites  using  different  approaches, 
where  little  effect  of  particle  clustering  was  observed  [2,55]. 


5.4.  Effect  of  spatial  homogeneity  on  plastic  behavior 

In  contrast  to  the  relatively  minor  effect  of  microstructural 
homogeneity  on  the  elastic  response  of  the  composite  mod¬ 
els,  there  is  a  much  stronger  influence  of  micro  structure  on  the 
plastic  behavior.  Fig.  7  shows  the  composite  yield  strength  mea- 


Also  shown  are  data  from  the  Eshelby  predictions  and  finite  element  results  for 
the  “homogeneous”  composite  models,  for  comparison. 
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sured  in  both  the  axial  and  transverse  directions  plotted  versus 
and  •  Since  we  are  restricting  our  modeling  to  the  onset 
of  plastic  flow  (i.e.  total  strains  <1%),  we  choose  as  the  com¬ 
posite  yield  strength  an  offset  yield  stress  corresponding  to  the 
flow  stress  at  an  overall  plastic  strain  of  2  x  10  4  (i.e.  cty0'02%). 
This  allows  a  direct  comparison  of  the  different  micro  structures 
beyond  the  micro-yielding  regime,  but  before  significant  strain¬ 
hardening  has  occurred  [56]. 

There  is  clearly  a  much  stronger  effect  of  microstructural 
heterogeneity  on  the  yielding  behavior  of  the  composites  than 
is  observed  for  the  elastic  behavior,  for  example  there  is  a 
difference  of  20-30  MPa  between  the  yield  strengths  of  the 
least  homogeneous  microstructure  and  the  most  homogeneous 
microstructure  measured  in  both  the  axial  and  transverse  direc- 
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tions.  Also,  there  is  slightly  less  scatter  in  the  yield  strength 
data  when  it  is  plotted  versus  the  transverse  homogeneous  length 
scale,  ,  compared  to  when  it  is  plotted  versus  the  axial  homo¬ 

geneous  length  scale,  .  This  may  suggest  that  this  particular 
microstructural  descriptor  is  slightly  more  relevant  than  its  coun¬ 
terpart  in  the  axial  direction.  This  observation  will  be  discussed 
further  in  the  context  of  strain  localization  in  the  inter-cluster 
regions  (see  Section  6). 

Additional  observations  regarding  the  plastic  behavior  of 
the  composites  are  presented  in  Fig.  8  which  shows  the  effect 
of  microstructural  homogeneity  on  the  initial  strain-hardening 
response  of  the  different  composite  microstructures.  Also  plot¬ 
ted  for  comparison  are  data  for  the  strain-hardening  rates  cal¬ 
culated  for  the  homogeneous  composite  models  and  predictions 
from  the  Eshelby  model.  The  strain-hardening  rate  is  evaluated 
at  overall  strains  of  between  0.005  and  0.006. 

From  the  figures  it  is  apparent  that  there  is  a  strong  effect 
of  microstructural  homogeneity  on  the  initial  strain-hardening 
behavior  of  the  composite  microstructures.  This  effect  is  clearly 
stronger  than  the  effect  on  yield  strength,  which  is  in  turn 
much  stronger  than  the  effect  of  microstructural  homogeneity 
on  elastic  modulus.  Furthermore,  there  is  a  stronger  effect  of 
microstructural  homogeneity  on  the  transverse  response  than 
the  axial  response.  For  example,  there  is  a  factor  of  around 
1 .4  between  the  initial  strain-hardening  rates  calculated  for  the 
axial  direction  in  the  most  homogeneous  and  least  homogeneous 
composite  microstructures.  This  factor  more  than  doubles  when 
comparing  the  same  microstructures  in  the  transverse  direc¬ 
tion,  suggesting  a  higher  sensitivity  to  microstructural  homo¬ 
geneity  in  this  loading  direction.  This  significant  knock-down 


Fig.  8.  Finite  element  results  for:  (a)  axial,  and  (b)  transverse  composite 
strain-hardening  rates,  plotted  vs.  the  microstructural  descriptors  and 
Microstructural  homogeneity  has  a  positive  effect  on  the  strain-hardening  rates 
evaluated  in  both  loading  directions,  with  a  stronger  effect  being  observed  for 
loading  in  the  transverse  direction  compared  with  the  axial  direction. 


Fig.  9.  Contour  plot  showing  regions  of  high  equivalent  plastic  strain  in 
PSR  =  8:1  composite  microstructure,  obtained  under  transverse  loading  at  a 
far-held  applied  stress  of  451  MPa  for  material  with  PSR  =8:1  under  trans¬ 
verse  loading.  The  maximum  (white)  contour  corresponds  to  =  0. 15,  loading 
direction  is  vertical. 
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effect  on  the  strain-hardening  ability  of  these  composites  could 
have  important  ramifications  in  terms  of  ultimate  strengths  and 
the  failure  mechanisms  operating  under  different  conditions  of 
homogeneity. 

6.  Discussion 

In  order  to  better  understand  the  role  of  microstructural  homo¬ 
geneity  on  determining  the  plastic  behavior  of  the  different 
composite  microstructures,  contour  plots  of  stress  and  strain 
field  quantities  were  generated  from  the  FEM  output  files.  For 
example,  Fig.  9  shows  an  elongated  band  of  high  equivalent 
plastic  strain  in  the  PSR  =  8:1  material  subject  to  transverse 
loading  (vertical  on  the  figure).  There  is  intense  strain  local¬ 
ization,  wherein  the  most  intense  plastic  straining  is  localized 
to  regions  of  lower-than- average  volume  fraction.  The  existence 
of  these  regions  is  a  direct  consequence  of  increasing  local  vari¬ 
ability  of  reinforcement  fraction  with  increasing  microstructural 
heterogeneity,  CV(Af).  For  a  constant  average  reinforcement 
fraction,  the  formation  of  higher- than- average  volume  fraction 
regions  (i.e.  clusters)  necessitates  the  associated  formation  of 
regions  of  low  volume  fraction,  i.e.  the  inter-cluster  regions.  The 
FEM  results  suggest  that  the  clusters  play  a  role  in  controlling 
deformation,  by  providing  local  constraint  to  the  plastically- 


deforming  regions.  For  example,  the  three  circled  regions  in 
Fig.  10(c)  are  regions  that  are  experiencing  enhanced  localized 
plastic  deformation,  due  to  the  interactions  between  the  adjacent 
high  volume  fraction  regions  (the  clusters)  and  the  applied  load, 
which  is  horizontal  in  the  figure. 

The  spatial  arrangement  of  the  clusters  with  respect  to  the 
loading  axis  is  therefore  important:  clusters  arranged  in  pairs 
parallel  to  the  applied  load  have  a  higher  propensity  for  strain 
localization  in  the  inter-cluster  region  than  those  oriented  in  pairs 
normal  to  the  applied  load.  Microstructural  anisotropy  (i.e.  dif¬ 
ferences  between  and  L^J)  can  therefore  play  an  important 
role  in  controlling  the  response  of  these  materials.  To  illustrate 
this,  Fig.  11  shows  the  data  for  initial  strain-hardening  rates 
measured  in  each  of  the  microstructures  re-plotted  versus  the 
homogenous  length  scale  measured  perpendicular  to  the  load¬ 
ing  axis.  In  the  plot,  both  the  axial  and  transverse  loading  data 
are  combined,  i.e.  for  axial  loading,  doq/dsi  is  plotted  versus 
and  for  transverse  loading,  d<72/d£2  is  plotted  versus  . 

The  figure  reveals  a  single  relationship  for  both  axial  and 
transverse  loading,  between  the  initial  strain-hardening  rate  and 
the  homogeneous  length  scale  measured  perpendicular  to  the 
loading  axis.  Increasing  the  microstructural  homogeneity  in  the 
direction  perpendicular  to  the  load  results  in  much  higher  strain¬ 
hardening  rates.  This  effect  can  be  explained  by  considering  the 


Fig.  10.  Indications  of  plastic  strain  localization  (circled)  in  PSR  =  8:1  microstructure:  (a)  binary  image  of  particles,  (b)  finite  element  mesh  of  same  region,  (c) 
equivalent  plastic  strain  contours.  Loading  is  applied  axially  (i.e.  horizontal)  with  a  far-field  stress  of  451  MPa. 
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Fig.  11.  Plot  showing  initial  strain-hardening  rates  d<xi/d£i  and  d(T2/de2  re_ 
plotted  vs.  the  homogeneous  length  scale  measured  transverse  to  each  loading 
axis  (i.e.  compare  with  Fig.  8(a)  and  (b)).  The  strain-hardening  data  all  collapses 
onto  a  single  line,  suggesting  that  there  is  a  general  relationship  between  the 
initial  strain-hardening  rate  and  the  spatial  homogeneity  of  the  microstructure, 
measured  perpendicular  to  the  applied  load. 


influence  of  microstructural  homogeneity  on  the  formation  of 
extended  regions  of  localized  plastic  strain,  analogous  to  the 
deformation  bands  which  are  familiar  in  homogeneous  materi¬ 
als  that  exhibit  very  low  work-hardening  rates  [11].  For  example, 
under  axial  loading,  it  is  the  spatial  homogeneity  in  the  trans¬ 
verse  direction  (i.e.  Lm)  which  is  the  relevant  microstructural 
descriptor.  As  Lh2  increases  (i.e.  it  becomes  less  homogeneous), 
the  probability  for  forming  extended  regions  of  low  reinforce¬ 
ment  volume  fraction  also  increases.  This,  in  turn,  facilitates 
strain  localization  in  these  regions  which  leads  to  a  lower  overall 
strain-hardening  response  in  the  composite.  Similar  arguments 
hold  for  loading  in  the  transverse  direction,  with  the  added  caveat 
that  in  these  composite  micro  structures,  the  spatial  homogeneity 
in  the  axial  direction  (i.e.  perpendicular  to  the  transverse  load¬ 
ing  axis)  is  typically  lower  than  in  the  transverse  direction  due 
to  microstructural  anisotropy  imparted  via  extrusion  process¬ 
ing;  Lhi  >Lh2-  Consequently,  strain-hardening  rates  measured 
in  the  transverse  direction  should  be  lower  than  those  measured 
in  the  axial  direction.  This  is  clearly  demonstrated  in  the  strain¬ 
hardening  data  presented  here. 

In  general,  there  is  a  significant  increase  in  the  amount  of 
strain  localization  when  transverse  loading  response  is  compared 
with  the  axial  loading  response.  In  order  to  show  quantitatively 
the  effect  of  microstructural  homogeneity  on  the  distribution  of 
stresses  and  strains  within  the  composites,  mean  equivalent  plas¬ 
tic  strains  were  computed.  Table  3  shows  the  mean  equivalent 
plastic  strains  generated  at  far-field  applied  stresses  of  45 1  MPa, 
for  both  axial  and  transverse  loading  directions.  This  particular 
stress  level  was  chosen  to  reflect  the  post-yielding  behaviors  of 
each  microstructure.  The  table  shows  that  for  an  applied  stresses 
of  45 1  MPa  in  the  axial  loading  direction,  there  is  an  increase  in 
plastic  strain  as  microstructural  homogeneity  decreases.  In  the 


Table  3 


Mean  values  of  equivalent  plastic  strain  generated  under  axial  and  transverse 
loading  with  a  far-field  applied  stress  of  45 1  MPa 


Equivalent  plastic  strain, 

£pi  (x  io-6) 

PSR  =  2:1 

PSR  =  3:1 

PSR  =  8:1 

Homogeneous 

Axial  loading 

207 

239 

437 

0a 

Transverse  loading 

301 

480 

2709 

0a 

Also  included  are  mean  equivalent  plastic  strains  for  a  “homogeneous”  compos¬ 
ite  with  an  overall  volume  fraction  of  0.275. 


Indicates  no  yielding. 


transverse  loading  direction,  this  effect  is  even  more  dramatic. 
The  rapidly  increasing  plastic  strain  is  due  to  the  strain  partition¬ 
ing  that  is  occurring  during  deformation.  For  example,  Fig.  12 
shows  that  there  is  a  dramatic  increase  in  the  amount  of  plastic 


(a)  element  volume  fraction  (-) 


(b)  element  volume  fraction  (-) 


Fig.  12.  Histograms  showing  the  distribution  of  equivalent  plastic  strain  vs. 
element  volume  fraction:  (a)  axial  loading,  (b)  transverse  loading.  Equivalent 
plastic  strains  were  generated  at  a  far-field  stress  of  45 1  MPa.  There  is  a  dramatic 
increase  in  plastic  strains  for  elements  with  lower  volume  fractions,  and  this  is 
further  enhanced  as  microstructural  homogeneity  decreases,  due  to  intense  strain 
partitioning. 
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strain  seen  in  elements  with  low  volume  fractions.  Furthermore, 
there  is  a  clear  effect  of  microstructure,  with  the  least  homo¬ 
geneous  material  exhibiting  the  highest  levels  of  plastic  strain 
at  the  lower  volume  fractions,  due  to  increased  levels  of  strain 
partitioning  in  these  materials.  The  effect  is  most  severe  in  the 
case  of  PSR  =  8:1,  loaded  in  the  transverse  orientation,  where 
ep?  ~  0-005  for  the  lowest  volume  fraction  elements,  whereas 

the  average  is  only  ~0.00089  (i.e.  an  increase  of  5.6x). 
At  the  other  end  of  the  scale,  the  “homogeneous”  composite 
microstructure  shows  zero  plastic  strain  at  this  level  of  applied 
stress,  i.e.  the  material  has  not  yet  yielded. 

If  these  strain  localization  effects  can  be  extrapolated  to  much 
higher  strain  levels,  the  relationships  between  spatial  homo¬ 
geneity  and  strain  localization  identified  by  the  present  model 
may  also  hold  implications  for  the  potential  failure  mechanisms 
which  can  operate  in  these  materials.  In  the  case  of  extremely 
low  spatial  homogeneity  (e.g.  high  PSR),  the  propensity  for 
strain  localization  is  high,  and  the  projected  strain-hardening 
rate  is  low,  therefore,  it  could  be  argued  that  a  strain-controlled 
failure  mechanism  such  as  void  nucleation  and  coalescence 
might  occur.  Conversely,  for  more  homogeneous  composite 
microstructures  wherein  strain  localization  is  suppressed,  a 
stress-controlled  mechanism  such  as  particle  cracking  and/or 
interfacial  debonding  might  operate.  In  the  former  case,  the 
strain-hardening  behavior  of  the  matrix  itself  would  play  an 
important  role  in  determining  the  overall  failure  strain  of  the 
composite.  Increasing  the  spatial  homogeneity  of  the  composite 
microstructure  is  could  therefore  be  considered  a  potential  way 
to  improve  the  mechanical  behavior  of  the  material,  by  allowing 
other  microstructural  enhancements  such  as  stronger  particles 
[57,58]  and/or  stronger  interfacial  bonding  [59,60]  to  dominate 
the  mechanical  response  of  the  composite,  rather  than  the  prop¬ 
erties  of  the  matrix. 

7.  Conclusions 

The  effect  of  spatial  homogeneity  on  the  elastic-plastic 
mechanical  response  of  extruded  Al-6061/SiC/27.5p  DRA  has 
been  predicted,  using  a  combined  approach  for  microstruc¬ 
tural  characterization  and  modeling.  Spatial  homogeneity  is 
quantified  using  an  extended,  anisotropic  Multi-Scalar  Analysis 
of  Area  Fractions  procedure,  resulting  in  two  microstructural 
descriptors:  the  homogeneous  length  scales,  Lhi  and  Lh2,  mea¬ 
sured  parallel  and  transverse  to  the  extrusion  axis,  respectively. 
A  new  model  for  predicting  the  effect  of  spatial  homogeneity 
on  mechanical  behavior  was  proposed  which  sought  to  incor¬ 
porate  a  mean-field  (Eshelby)  micromechanical  model  within 
a  simple  finite  element  model,  thereby  dramatically  increas¬ 
ing  the  model’s  extent  without  increasing  its  complexity.  This 
model  enables  the  investigation  of  deformation  mechanisms 
occurring  over  length  scales  encompassing  many  particle  diame¬ 
ters,  such  as  plastic  strain  localization  and  deformation  banding, 
without  requiring  long  run-times,  by  using  a  volume- averaged, 
discretized  approach.  A  key  assumption  in  this  approach  is  that 
below  the  length  scale  of  the  discretization,  the  microstructure 
can  be  modeled  by  a  uniform  array  of  reinforcement  parti¬ 


cles.  The  effect  of  the  level  of  discretization  on  predictions  of 
microstructure-property  relationships  are  not  addressed  directly 
in  the  present  work,  however,  detailed  comparisons  of  the  present 
model  with  discretely-modeled  matrix-particle  geometries  will 
form  the  basis  of  a  subsequent  publication,  to  encompass  phe¬ 
nomena  such  as  localization  of  plastic  deformation  between 
particles,  and  subsequent  matrix-particle  debonding.  Neverthe¬ 
less,  empirical  models  were  generated  which  describe  the  effect 
of  spatial  homogeneity  on  the  elastic  and  plastic  behavior  of  the 
composite  microstructures.  The  results  confirmed  that  spatial 
homogeneity  plays  only  a  minor  role  in  elastic  deformation,  but 
has  a  larger  effect  on  the  yield  strength.  By  far  the  strongest  effect 
is  on  the  strain-hardening  rate,  where  the  formation  of  deforma¬ 
tion  bands  is  promoted  by  the  more  heterogeneous,  anisotropic 
arrangements  of  the  SiC  particles.  There  is  a  strong  directionality 
to  this  effect,  with  the  relevant  (i.e.  controlling)  microstructural 
descriptor  identified  as  the  homogeneous  length  scale  measured 
perpendicular  to  the  loading  direction.  A  universal  relationship 
is  proposed  between  the  homogeneous  length  scale  measured 
perpendicular  to  the  loading  direction  and  the  strain-hardening 
rate.  It  is  postulated  that  the  key  mechanism  controlling  strain¬ 
hardening  behavior  is  the  interruption  of  the  deformation  bands 
which  form  in  elongated  low  volume  fraction  regions  in-between 
the  SiC  clusters.  Results  from  the  FEM  analysis  suggest  that  as 
spatial  homogeneity  decreases,  there  is  a  steady  increase  in  the 
plastic  strain  in  these  regions;  equivalent  plastic  strains  can  reach 
5-6  times  the  average  value  for  the  most  clustered  microstruc¬ 
tures,  at  the  chosen  level  of  discretization.  This  effect  is  more 
noticeable  for  the  transverse  loading  case,  due  to  the  anisotropy 
of  the  extruded  microstructures. 
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